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ABSTRACT 


We report the first detection of a time-dependent, dust-scattered X-ray halo 
around a gamma-ray burst. GRB 031203 was observed by XMM-Newton starting 
six hours after the burst. The halo appeared as concentric ring-like structures 
centered on the GRB location. The radii of these structures increased with time 
as t 1 / 2 , consistent with small-angle X-ray scattering caused by a large column of 
dust along the line of sight to a cosmologically distant GRB. The rings are due 
to dust concentrated in two distinct slabs in the Galaxy located at distances of 
880 and 1390 pc, consistent with known Galactic features. The halo brightness 
implies an initial soft X-ray pulse consistent with the observed GRB. 
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1. Introduction 

It has long been realized that the small-angle scattering of X-rays by dust grains can 
result in a detectable X-ray “halo” around a distant X-ray source (Overbeck 1965), with a 
radial intensity distribution that depends on the dust properties and location. This phe- 
nomenon was first detected by Rolf (1983) using data from Einstein and confirmed by later 
observations (e.g. Mauche & Gorenstein 1986; Predehl et al. 1991). 

In the case of Gamma- Ray Bursts (GRBs) , the transient nature of the burst combined 
with the large initial X-ray flux means X-ray scattering may be visible for a short period of 
time with a relatively high surface brightness. As the transient flux passes through a dust 
slab between the observer and the GRB, the scattering process introduces a time-delay in 
which the X-rays at larger angles to the line of sight arrive increasingly delayed with respect 
to the non-scattered X-rays. Thus, the dust allows us to view the GRB X-ray flux at earlier 
times. The time-dependent X-ray halo around the GRB can, in principle, be used to provide 
detailed information on the location, spatial distribution and properties of the dust and 
the distance to and brightness of the GRB (e.g. Triimper & Schofelder 1973; Klose 1994; 
Miralda-Escude 1999; Draine 2003). Images of transient X-ray sources seen in scattered light 
are analogous to “light echos” seen around some supernovae (e.g. SN 1987A, Xu, Crotts & 
Kunkel 1995). 

Scattered X-rays observed at an angle 9 from the line-of-sight to a GRB arrive after the 
direct emission with a time-delay t s given by: 



where z s and D s are the redshift and angular diameter distance of the scattering dust, 
D g is the angular diameter distance of the GRB and D gs is the angular diameter distance 
from the dust to the GRB (Miralda-Escude 1999). 

Here we discuss the case of GRB 031203 observed by XMM- Newton on two occasions 
shortly after the burst. The first observation revealed the first dust-scattered X-ray halo 
detected around a GRB. 


2. Observations and data analysis 

GRB 031203 was detected by the IBIS instrument on Integral on 2003 December 3 at 
22:01:28 UT as a single peaked burst with a duration of 30 s and peak flux of 1.3 x 10 -7 
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erg s -1 cm -2 in the 20-200 keV band (Gotz et al. 2003; Mereghetti &: Gotz 2003). A 58 ks 
XMM-Newton observation of the field began at 2003 December 4, 04:09:29 UT. Several X-ray 
sources were detected within the 2.5' radius Integral error-circle (Santos-Lleo & Calderon 
2003), the brightest of which (Si) appeared to fade through the observation (Rodriguez- 
Pascual et al. 2003) and was interpreted as the X-ray afterglow of GRB 031203. We refined 
the astrometric position of the afterglow by performing a cross-correlation with the USNO- 
A2 catalog, the improved Si position was (J2000) RA=08 h 02 m 30.19 s , Dec= —39° 51' 04.0" 
(1<7 error 0.7", Tedds et al. 2003). 

The position in the Galactic plane (l = 255°, b = —4.6°) results in high optical extinc- 
tion (E(B — V) = 1.0, Schlegel, Finkbeiner & Davis 1998). Optical/IR observations failed 
to locate an afterglow. Radio observations did locate a transient source location consistent 
with the X-ray afterglow (Soderberg, Kulkarni & Frail 2003), and optical imaging and spec- 
troscopy revealed a candidate host galaxy at 2 = 0.105 which shows star-forming features 
typical of GRB hosts (Prochaska et al. 2003a, b). Our refined location for the X-ray afterglow 
is within 0.5" of the optical galaxy and consistent with the position of the radio source. 

During the first XMM-Newton observation the GRB was observed on axis with all EPIC 
instruments in full-frame modes (Striider et al. 2001; Turner et al. 2001). A second 54 ks 
XMM-Newton observation was obtained on 2003 December 6 in which the GRB was « 6' 
off axis. The observational details and the GRB spectrum and light-curve are discussed in 
detail 3 * 5 in Watson et al. (2004). 

The extraction of science products followed standard procedures using the XMM-Newton 
Science Analysis System (SAS v5 . 4 . 1). The EPIC data were processed using the SAS chains 
to produce calibrated event lists and remove background flares. GRB afterglow source counts 
were extracted from circular regions (radius 34") and background counts were estimated from 
a large off-axis region free from obvious point sources and the halo. The GRB afterglow 
fades as (t — t Q )~ 0A , where t 0 is the time of the GRB, which is unusually slow for GRB X-ray 
afterglows. 


3. The X-ray Halo 

The X-ray image from the first XMM-Newton observation reveals an extended, circular 

halo concentric with the GRB 031203 afterglow (see Fig. 1). This halo was seen in all three 
cameras of the EPIC instrument, is unique to this observation, and is not due to scattered 
optical or X-ray light within the instrument. Summing the entire observation the halo had 
the form of a virtually complete ring (Vaughan et al. 2003), but splitting the observation 
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into contiguous time intervals revealed distinct rings which increased in radius through the 
observation. 

In order to better quantify this expansion, EPIC MOS images were produced for ten 
contiguous time intervals of 5780 s duration, spanning the entire first observation, and binned 
to 4" pixels. We concentrate on the MOS data as the halo lies over several chip gaps in the 
EPIC pn images, but the same halo was seen in both. The radial profiles of the images about 
the afterglow were calculated removing the three brightest point sources within the region 
of interest (SI, S2 and S3 in Fig. 1) by ignoring all counts within 40" from the centroid of 
each source. The radial profiles were then background subtracted using the mean level from 
the 350 - 400" annulus, which is well outside the detected halo. The resulting count profiles 
as a function of radius are shown in Fig. 2. The strongest peak moved outward from ~ 120" 
to ~ 220" during the observation, while a second peak moved from ~ 160 to ~ 270". 

The radial positions of the two rings were measured from the local maxima of the two 
peaks in the radial distributions. Fig. 3 shows the change in radii as a function of time, 
and both were well fit by a simple power-law (9 oc (t — to)“)- Fixing to — 0, the indices 
of the power-law expansions were measured to be a = 0.53 ± 0.04 and 0.51 ± 0.05 for the 
inner and outer rings respectively (errors are 90% confidence intervals), consistent with the 
9 oc (t — 1 0 ) 1/2 expansion predicted for a scattered halo. As a consistency check we allowed 
t 0 to be a free parameter but fixed a — 0.5 and refitted to derive limits of t 0 = 2794^3^8 s 
and 2005^2867 s f° r the inner and outer ring respectively, both consistent with t 0 = 0. The 
second XMM-Newton observation obtained 3 days after the burst showed no evidence for 
any ringed structure. 

The distance to the GRB is large (adopting the redshift of the putative host galaxy), 
the time-delay is short and 9 is relatively large, thus to a good approximation D g — D gs and 
(1 + z s ) = 1 in equation (1), hence D s = 2 ct s /9 2 . Assuming 9 oc (t — to) 1 / 2 and t 0 = 0, the 
data in Fig. 3 imply that the scattering dust slabs are located at distances from the observer 
of £>i = 1388 ± 32 pc and D 2 = 882 ± 20 pc, corresponding to the inner and outer ring 
respectively. At these distances, the ring diameters probe size scales of 2-3 pc in the dust 
scattering medium. 

The halo spectrum was extracted from an annulus with radii 110 — 220" using only the 
first half of the observation; the decreased surface brightness of the halo at later times makes 
spectral extraction unreliable. A background spectrum was extracted from a large off-axis 
region (avoiding the halo and other point sources). Fig. 4 shows the halo spectrum compared 
to the afterglow spectrum. The spectrum of the halo is much steeper than that of the 
afterglow. Fitting with a simple absorbed power-law yielded photon indices of T = 1.98±0.05 
and T = 3.03 ± 0.14 for the afterglow and halo, respectively. The best-fitting absorbing 
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column was N-r = 8.8 ± 0.5 x 10 21 cm -2 when fixed to be the same for both afterglow and 
halo spectra. This is larger than the Galactic value of Nr = 6.1 x 10 21 cm -2 (Dickey & 
Lockman 1990). The specific fluxes at t = 36650 s after the burst at 1 keV are 1.5 +q 2 x 10~ 4 
photons cm -2 s -1 keV -1 for the afterglow and 6.0 ± 0.4 x 10 -4 photons cm -2 s -1 keV -1 for 
the halo. 

To model the dust scattering medium, we have simultaneously fitted the ten radial pro- 
files shown in Fig. 2 using the differential and total scattering cross sections of the Rayleigh- 
Gans approximation (Mauche and Gorenstein 1986). This is a valid approximation if the 
grain radius in pm is much less than the X-ray photon energy in keV (Smith & Dwek 1998). 
We assume that the dust is confined to two slabs at distances D\ and D 2 with thickness AD. 
The halo seen is created by the combination of the soft X-ray pulse from the GRB and the 
subsequent afterglow. The observed angular width of the rings results from a combination 
of the duration of the soft X-ray pulse, At, AD and the point spread function (PSF) of the 
telescope. The width of the PSF corresponds to At ~ 1000 s or AD ~ 100 pc. Given that 
the GRB is very short, ~ 30 s (Mereghetti & Gotz 2003), we assume that At <C 1000 s and 
hence, allowing for the PSF, the observed broadening of the rings, 20", gives us an estimate 
of the slab thickness, AD = 130 ± 50 pc. 

The angular and temporal distribution was well fitted using grains with radii in the 
range 0.15 < a < 0.25 pm. Smaller grains scatter to much larger angles (not visible because 
the time elapsed since the burst is too small) and the upper limit to the grain size distribution 
is expected to be « 0.3 pm (Predehl et al. 1991). Using this range of grain sizes, we require 
the initial GRB pulse and afterglow to have a photon index of T ~ 2.0 (as observed for the 
afterglow from 6 h) to match the observed photon index of the halo (F = 3.0). 

The brightness of the halo depends on the product of the differential scattering cross 
sections of the grains, the column densities of the grains in the slabs and the integrated flux 
of the X-ray pulse and afterglow. Interstellar extinction maps show structure in the GRB 
direction, with a considerable increase to Ay « 2 mag at a distance of ~ 1.3 kpc (Neckel & 
Klare 1980). This distance is consistent with the location of the more distant dust slab we 
detect. Adopting a mean grain size of a = 0.2 pm, an Ay = 2 mag corresponds to a grain 
column density of 1.5 x 10 8 cm -2 (Mauche and Gorenstein 1986). The best fit to the ring 
structure (x 2 /d.o.f. = 836/489) results from dividing this column density between the two 
slabs, TVi ~ 1 x 10 s and N 2 ~ 0.5 x 10 8 cm -2 (i.e. most of the dust is in the more distant 
slab). Using these values and a = 0.20 ± 0.05, the time-integrated flux required for the soft 
X-ray pulse is 1600 ± 800 photon cm -2 keV -1 at 1 keV. Assuming a 30 s rectangular pulse 
this is equivalent to ~ 5 Crabs; using T = 2 the predicted GRB flux in the 20 — 200 keV 
band is 2.4 ±1.2 photon cm -2 s -1 . Given the assumptions (e.g. spectral shape) this is in 
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reasonable agreement with the value reported for GRB 031203 (peak of ~ 1.2 photon cm -2 
s _1 , Gotz et al. 2003; Mereghetti & Gotz 2003). The X-ray burst is very bright compared to 
the afterglow (see Watson et al. 2004). Extrapolating the decay of the afterglow (section 2) 
to the time of the GRB accounts for only ~ 2 per cent of the halo emission. 

The line of sight to GRB 031203 passes close to the center of the Gum Nebula which 
appears as a 28° diameter sphere in Ho: centered at l = 258°, b = —2° (Chanot & Sivan 
1983). The Gum Nebula is likely to be a supershell created by repeated SNe explosions 
within it. Reynoso & Dubner (1997) detected a neutral gas disk (H i) associated with the 
Gum Nebula, located ~ 500 ± 100 pc from us, with a radius of ~ 150 pc. The dust slab 
nearest to us may be associated with the rear face of the Gum Nebula which is known to 
have a concentration of molecular clouds (Woermann, Gaylard & Otrupcek 2001). 


This paper is based on observations obtained with XMM-Newton, an ESA science mis- 
sion with instruments and contributions directly funded by ESA Member States and the USA 
(NASA). SV and AJL acknowledge support from PPARC, UK. We acknowledge the ben- 
efits of collaboration within the Research and Training Network “Gamma-Ray Bursts: An 
Enigma and a Tool”, funded by the European Union under contract HPRN-CT-2002-00294. 
We thank an anonymous referee for a constructive report. 
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Fig. 1. — Combined MOS images of GRB 031203 covering the 0.7 — 2.5 keV range, spanning 
10' on a side. The data were divided into 10 contiguous time intervals lasting 5780 s. The 
images were obtained from the first six time slices are shown (smoothed using a 6 arcsec 
Gaussian kernel). The three brightest point sources (Si, S2, S3) are marked. 

Fig. 2. — Radial profile of counts about the GRB 031203 afterglow from ten contiguous 
time intervals (derived from the combined MOS 0.7 — 2.5 keV images). The afterglow itself 
(Si) and the next two brightest points sources (S2, S3) were removed prior to calculating 
the profiles (the positions of S2 and S3 are marked in the upper panels). The two peaks 
corresponding to the two expanding rings are obvious. The dotted lines denote the best- 
fitting 9 oc (t — to) 1 / 2 functions (see Fig. 3). 

Fig. 3. — Expansion of the two rings around GRB 031203 with time. The radius of each 
ring was measured from the local maxima of the radial profiles (Fig. 2). In both cases the 
expansion was well-fitted with a functional form 9 oc (t — to) 1 / 2 . The errors were assumed to 
be the width of the bin size used for the radial profile measurement. 

Fig. 4. — Absorbed power-law model fits to the EPIC pn spectra of the GRB 031203 afterglow 
(circles, thick line) integrated over the entire observation and the dust halo (squares, thin 
line) integrated over the first half (27895 s) of the observation (during which the surface 
brightness was high enough for a reasonable spectral extraction). This clearly demonstrates 
the spectrum of the halo is significantly steeper than the afterglow. 
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